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Abstract 

In this article we present triple differential cross sections for equal energy sharing 
kinematics for double photoionization of the helium atom at 20 and 40eV above 
threshold in the framework of the hyper spherical partial wave theory. This sup- 
plements our earlier work [1] in which we were successful in showing fully, gauge 
independence of the results in our formalism. Also in this work we treat cases in 
which the Stokes parameter Si < 1 so that partial polarization of the photon source 
is also taken into account. Agreement in shape with the Convergent Close Coupling 
[14,17] calculation and the experiments appears to be excellent. 
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1 Introduction 



Double photo ionization (DPI) of the helium is one of those fundamental 
atomic processes which has received wide attention from both theorists and 
experimentalists. Helium is one of the most abundant elements in stellar atmo- 
spheres and hence analysis of Helium spectra is of considerable importance in 
astrophysics. From the theoretical aspect, the absorption of a photon and sub- 
sequent ejection of both electrons from the Helium atom into the continuum, 
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leads to a pure Coulomb three body problem which, at low incident energies, 
involve complex short range interactions and strong electron correlations that 
are difficult to treat theoretically and hence present a challenge to theorists. 

Detailed discussion about existing theories of helium double photo ionization 
is presented in our earlier work [1]. It would still be worthwhile to mention 
the convergent close coupling (CCC) theory of Bray et al [2], 2SC theory of 
Pont and Shakeshaft [3] and the 3C theory of Maulbetsch and Briggs [4] since 
all of them have been very widely used to study the double photo ionization 
process. It would also be very appropriate to single out and mention the hyper- 
spherical 7^-matrix with semi classical outgoing waves (H7?.M-S0W) theory of 
Malegat et al [5,6] for two reasons. First, the method has been quite successful 
in describing the total ionization cross section (TICS), single and triple differ- 
ential cross sections (SDCS and TDCS). Secondly the method bears similarity 
to our calculation in the use of hyperspherical coordinates, though the tech- 
niques used in extracting scattering information is very different from ours. 
The H7?.M-S0W method uses 7^-matrix [7] method to obtain solution of the 
radial wave equation in the inner region bounded by a hypersphere at Rq {R 
being the hyperradius) . Beyond Rq the solution of the radial equations are 
propagated semiclasically (see section IV of [6]). In contrast our method is 
fully quantum mechanical and free from any approximations, though perhaps 
it is more demanding on the computational resources. 

The helium double photo ionization is similar to the electron impact ionization 
of atomic hydrogen. However, due to the ^S^ symmetry of the initial state, the 
final state can be only of symmetry. This leads to considerable simplifica- 
tion in computations in contrast to e-H ionization, in which there are many 
more contributing states in the final channel. Also there are selection rules 
that are valid only for double photo ionization [4] (for a recent review see [8]). 

As is well known, the DPI triple differential cross-section can be calculated in 
three gauges, namely length, velocity and acceleration gauges. In principle the 
TDCS must be independent of the choice of gauge. However, it has been shown 
by Lucy et al [9] that the TDCS results are notoriously gauge dependent un- 
less one uses very accurate wave functions for both initial and final channels. 
In our previous work [1] we had been able to show that our results are largely 
gauge independent and length and velocity gauge results obtained were nearly 
identical. We now have the TDCS results in the acceleration gauge also, and 
we will make a brief comparison between the TDCS in the three gauges below 
before we study other cases. Also in our work [1] we obtained TDCS at 20eV 
incident photon energy and our results were for incident photon source lin- 
early polarized with Stokes parameter 5*1 = 1. Having obtained an essentially 
gauge independent formalism, in this work we present double photo ioniza- 
tion TDCS at 20eV and 40eV above the helium double ionization threshold 
79eV. Also we choose situations in which the incident photon beam is partially 
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(linear) polarized, the degree of polarization being given by the Stokes 
parameter Si (for a definition of Stokes parameters see [8]. Our corre- 
sponds to Siin given in section 3.2 of Ref. [8] in which the x and y axes are 
respectively along the major and minor axes of the polarization ellipse) . 

The double ionization TDCS is obtained from the transition matrix elements 
given by 

T^, = (1) 



where $i(ri, r2) is the helium ground state, fi and r2 being the coordinates of 
the two electrons with respect to the helium nucleus (assumed to be at rest), 
is the final channel continuum state and is the interaction term given 

by 

Vi = e- D. (2) 



— * — * — * — * 

In (2) D is the dipole operator given hy D — Vi + V2 (velocity form), 

— * — * - — » — » 

D — uji{fi + r2) (length form) or D — — ^(Vi + V2)V (acceleration form), 
where 



ri r2 |fi - f2 



is the full three body interaction potential and e is the photon polarization 
vector. For an arbitrary degree of (linear) polarization of the incident photon 
beam characterized by the Stokes parameter Si, the TDCS is given by (see 
[19] for example and references therein) 



where ax and cxy refer to TDCS calculated with polarization vectors along the 
X-axis and y-axis respectively. 

A schematic diagram of the scattering geometry is shown in figure 1. The 
photon emerges from the bottom of the scattering plane along the z-axis. The 
two photo-electrons are ejected after collision in the direction 9i and 6*2 in the 
scattering plane. 



3 




Fig. 1. Schematic diagram of the scattering geometry. 

2 Theory 



The hyperspherical partial wave theory is discussed in considerable detail in 
our earlier work [1]. Therefore we give here only the essentials. 

For the Helium ground state we take a highly correlated 20 term Hylleraas type 
wave function given by Hart and Herzberg [10]. To calculate the final channel 

continuum state \E'j ■* we use hyperspherical coordinates R = \Jr'l + r^, a = 
arctan(r2/ri), fi = (6'i,0i), f2 = (6*2,02) and uj = (a,fi,f2). Also we set 

P = \/Pi+P2, «o = arctan(p2M), Pi = {9p^,(l)pJ, p2 = (Op^Api) and ujq = 
{ao,pi,p2) where and pi (i = 1, 2) are the coordinates and momenta of 
the i^^ particle, ■* is then expanded in hyperspherical harmonics [11,12] 
that are functions of the five angular variables and £i, £2, L, M, which are 
respectively the angular momenta of the two electrons, the order of the Jacobi 
polynomial in the hyperspherical harmonics, the total angular momentum and 
its projection. For a given symmetry s we decompose the final state as 

^^fs\R.^) = \I^T.^<Pli-) (4) 
V vr ^ p2 



where A is the composite index (£1, £2, n, L, M) or 2n + ii + £2 depending 
on the context and p = PR and (pxiuj) are are orthogonal functions that are 
a product of a Jacobi polynomial P^^'^(a) and a coupled angular momentum 
eigenfunction yl'J^ {il 1,^12) ■ 

then satisfy the infinite set of coupled differential equations 
2 , 1^. 2 at 



^ I^A (i^A + 1 

1 o ~r J- 



^ ,n(pn^^FM-0. (5) 



Here are the matrix elements of the full three-body interaction potential 
and i/^ = A + |. Since the final channel state must have the symmetry s 
is fixed. The contributing radial waves then correspond to L = 1 and (odd) 
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parity vr = — 1 so that writing = {^1,^2,^) and — 
the relevant set of radial waves become 



/jv the equations for 




/iv + E 



N' 



Pp 



0. 



(6) 



For actual computations we truncate the set of equations to some maximum 
value A^mx of N . These N„i,j. equations in N^^ variables are solved from origin 
to infinity. Construction of the radial wave solution is presented in our earlier 
works [1,13] with considerable rigor. So we omit the details in this work. 
Knowing the radial wave solution, the final channel state can be found from 
(4) and the transition matrix elements from (1). The photoionization TDCS 
can then be obtained using 



dQ.1dQ.2dE1 UJi 



3 Results and discussion 

Convergence in our cross sections depend on two parameters, namely the num- 
ber of coupled channels included in the computations and the asymptotic range 
parameter R^o where the asymptotic and the interior solutions are matched 
(see [1] and [13] for details). Numerical investigations show that convergence 
in TDCS is obtained at Roo = 30000 a.u. with 90 coupled channels. On the 
contrary oscillations in the SDCS persist even on increasing the asymptotic 
range parameter R^ substantially beyond 30000 a.u. This is in contrast to 
the HT^M-SOW method where oscillations in the SDCS die out on extending 
the asymptotic range. However, our computed SDCS at E/2 {E being the 
excess energy) remains more or less constant with the increase of convergence 
parameters. This strongly suggests that we employ a scaling technique similar 
to that used in the CCC theory (see [17]). Whenever the true SDCS at E/2 is 
known (from experiment or from some other theory, for example [15,16]) we 
normalize our computed SDCS at E/2 to the true SDCS by multiplying with 
a factor. Subsequently, we use this factor to scale our TDCS. Thus the TDCS 
at 20 eV incident energy have been scaled with a factor 0.8 and at 40 eV by 
a factor 0.6. In this way we are able to obtain absolute TDCS. 

3.1 Gauge independence 

As remarked earlier, in our previous work [1] we had been able to show that our 
work is largely gauge independent. There, independence of the TDCS results 



d^a 



27r2 



apiP2 




(7) 
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with respect to the length and velocity gauges had been presented. We now 
present the TDCS results in the acceleration gauge also. These are displayed 
in figure 2. Except near the peaks, the length and velocity gauge results are 
identical. The maximal departure of the length gauge results occur near the 
peaks, and are respectively 20% at 9i = 0", 14% at 9i = 30°, 14% at Oi = 60° 
and 7% at 9i = 90° from the velocity gauge. The acceleration gauge results 
are indistinguishable from that of the velocity gauge. The gauge independence 
in our results is a strong signature that our final channel wave function has 
the correct asymptotic and short range behaviour. Since the velocity and ac- 
celeration gauge results are identical, in subsequent calculations we compute 
TDCS in the acceleration gauge only, as this makes our computations much 
simpler for values of the Stokes parameter S*! < 1. 
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Fig. 2. TDCS for equal energy sharing double photo ionization of the helium atom at 
99eV incident photon energy in units of 10~^ a.u. and (a) 9i = 0°, (b) 6i = 30°, (c) 
9i = 60°, (d) 9i = 90°. Results are for length gauge (dashed curve), velocity gauge 
(continuous curve) and acceleration gauge (dotted curve). The Stokes parameter 
Si = 1.0. 
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3.2 Results for 20eV excess energy 





Fig. 3. TDCS for equal energy sharing double photo ionization of the helium atom 
at 99eV incident energy in units of 10^^ a.u. for (a) 6i = 180°, Si = 0.57, (b) 
$1 = -76°, Si = 0.56, (c) ei = -91°, = 0.53. Theory: continuous curve - present 
results; dashed curve - CCC results from [14]. In (b) and (c) the CCC results have 
been scaled with a factor 0.52. Experiment: Filled circles with errorbars are from 
[8] normalized suitably with the present results in each figure. 

In figure 3 we compare our results with experimental data from Briggs and 
Schmidt [8] and the results of the CCC theory [14,17] which are on an absolute 
scale. Since the experimental data are not absolute we have normalized the 
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data in each set to our computed results. Also the CCC results in figure 
3 (b) and (c) have been scaled with a factor 0.52. The agreement of the 
present results with the experiment and the CCC results appears to be good 
everywhere in shape, except for secondary peaks near 62 — ±150° at 6\ — 180° 
(figure 3 (a)). To see whether this additional structure in figure 3 (a) remains, 
we may need to perform a much larger calculation, that involves the inclusion 
of many more partial waves in our final state wave function. 
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Fig. 4. TDCS for equal energy sharing double photo ionization of the helium atom 
at 99eV incident energy in units of 10~^ a.u. Theory: (a) di = 98", (b) 9i = 115", 
(c) 61 = 132" The Stokes parameter = 0.67. Experiment: Filled circles with 
errorbars are from [18] normalized suitably with the present results by multiplying 
with a single factor. 

In figure 4 we compare our results with the experimental data of Weightman 
et al [18]. The experimental data is not absolute and have been normalized to 
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the present results by multiplying all data with a single factor. The agreement 
is excellent everywhere. 




Fig. 5. TDCS for equal energy sharing double photo ionization of the helium atom 
at 119eV incident energy in units of 10~^. Theory: (a) 9i = 95°, (b) Oi = 105°, 
(c) 6*1 = 115°, (d) 6*1 = 125°, (e) Oi = 135°. The Stokes parameter Si = 0.8. 
Experiment: Filled circles with errorbars are from [19] normalized suitably with 
the present results by multiplying with a single factor. The CCC results have been 
scaled by 0.4 for all figures. 



There is the slight hint of a peak in figure 4 (a) near 02 — 180° and its 

magnitude appears to increase in figures 4 (b) and (c). Again, confirmation 
of this or otherwise have to be deferred due to present limitations in our 
computational resources. The CCC results for this set of data are not available 
and hence comparisons are not made. 
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3.3 Results for 40eV excess energy 

Figure 5 shows our results for excess energy photon energy £^ = 40 eV and 
Stokes parameter 5*1 = 0.8. Shown also are the experimental results from 
Cvejanovic et al [19] and the CCC results presented in Ref. [19]. The exper- 
imental points have been scaled to our theory by multiplying with 0.05 in 
all cases. The absolute CCC results have also been scaled by the same factor 
0.4. Agreement in shape with the experimental results appear to be excellent 
everywhere, except in figures 5 (d) and (e) where there are slight departures 
around 62 = 150°. The CCC results are also almost identical with respect to 
the shapes. 



4 Conclusions 

In this work, we have presented results for equal energy sharing double photo 
ionization of the helium atom at 20 eV and 40 eV excess energy. Gauge inde- 
pendence of our TDCS results are shown. Cases in which the incident photon 
beam is partially polarized are considered. Comparisons are made with the 
experiments and the CCC theory wherever available and the results are seen 
to be consistent in shape. In the absence of absolute TDCS measurements for 
the chosen kinematics, it is difficult to say anything about the correctness of 
the magnitude of the various results. We also mention that our results pre- 
sented in this work have converged approximately. For getting fully converged 
results, more computational facilities may be necessary. All the computations 
reported here were done on desktop computers with Pentium IV class CPU 
and 512M core memory. 

In a future work, we propose to deal with unequal energy sharing kinematics. 
However, due to computational limitations we cannot reproduce results for 
extremely asymmetric energy sharing at this moment. As noted in our work 
[13] high Rydberg states tend to interfere with our continuum state giving 
undesirable results in such cases. To cope with these situations, considerably 
more computational resources may be necessary. 
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